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ABSTRACT 

It has been claimed that the observed magnitude of the vacuum energy density is 
consistent with the distribution predicted in anthropic models, in which an ensemble 
of universes is assumed. This calculation is revisited, without making the assumption 
that the CMB temperature is known, and considering in detail the possibility of a 
recollapsing universe. New accurate approximations for the growth of perturbations 
and the mass function of dark haloes are presented. Structure forms readily in the 
recollapsing phase of a model with negative A, so collapse fraction alone cannot forbid 
A from being large and negative. A negative A is disfavoured only if we assume that 
formation of observers can be neglected once the recollapsing universe has heated to 
r > 8 K. For the case of positive A, however, the current universe does occupy a 
extremely typical position compared to the predicted distribution on the A — T plane. 
Contrasting conclusions can be reached if anthropic arguments are applied to the 
curvature of the universe, and we discuss the falsifiability of this mode of anthropic 
reasoning. 
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1 INTRODUCTION 

There is now ahnost unanimous agreement that the universe 
contains a component that strongly resembles Einstein's cos- 
mological constant. This conclusion is cross-buttressed by a 
variety of lines of evidence (e.g. Spergel et al. 2006), and 
seems unlikely to change. Physically, this means that the 
vacuum appears to possess a homogeneous energy density, 
with an associated negative pressure. 

The existence of a non-zero vacuum density raises two 
problems: (1) the scale problem, and (2) the why-now prob- 
lem. The first of these concerns the energy scale correspond- 
ing to the vacuum density. If we adopt the values f2„ = 0.75 
and h — 0.73 for the key cosmological parameters, then 

p„ = 7.51x10-^^ kg m-3 = ^ (^)', (1) 

where Ev — 2.39 meV is known to a tolerance of about 
1%. The vacuum density should receive contributions of 
this form from the zero-point fluctuations of all quantum 
fields, and one would expect a net value for Ev of order 
the scale at which new physics truncates the contributions 
of high-energy virtual particles: anything from 100 GeV to 
10^^ GeV. The why-now problem further asks why we are 
observing the universe at almost exactly the unique time 
when this strangely small vacuum density first comes to 
dominate the cosmic density. 

Taking the second problem first, a question that in- 
volves the existence of observers must necessarily have an 
answer in which observers play a role. Therefore, a solu- 
tion to the why-now problem requires anthropic reasoning. 
In defining exactly what this term means, it is convenient 



to draw a distinction between what might be called one- 
universe and many-universe anthropic arguments. The first 
of these should not be controversial, since it states that ob- 
servers are likely to exist at special times in this universe: 
complex structures cannot begin to form until tempera- 
tures reach T <^ 1000 K, and the formation of structure 
will switch off in the near future when conditions become 
heavily vacuum dominated and the universe enters a phase 
of exponential expansion. Because the densities in baryons 
and dark matter are very roughly comparable, and because 
the photon-to-baryon ratio is roughly 137^mp/me, this win- 
dow for structure formation also in practice opens soon af- 
ter matter-radiation equality (e.g. p97 of Peacock 1999), al- 
though there is no suggestion that this is anything other 
than a genuine coincidence. If we further assume that com- 
plex life requires metals, it is no surprise that the universe 
now has an age comparable to a typical stellar lifetime, and 
this is almost certainly the explanation for Dirac's large- 
number coincidence (e.g. Carter 1974). 

The single-universe expectation that we should live rel- 
atively soon after matter-radiation equality may be all that 
is required of anthropic reasoning. The quintessence pro- 
gramme aims to find a dynamical origin for the vacuum 
energy; the hope is that the change in cosmological expan- 
sion history at matter-radiation equality may prompt the 
effective vacuum density generated by some scalar field to 
change from a sub-dominant contribution at early times, to 
something that resembles A by the present. This would be 
a satisfying solution to the why-now problem, but it is not 
clear that the mechanism can be made to work. For a sim- 
ple scalar field with an arbitrary potential, it seems that an 
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energy scale of ~ 1 meV needs to be introduced into the 
potential by hand, in order to prevent the quintessence den- 
sity from tracking the overall mass density at all times (e.g. 
Liddle & Scherrer 1999). This simply swaps one unexplained 
coincidence for another. Greater success is achieved with the 
more radical option of fc-essence, in which the field is given 
a non-canonical kinetic term. This readily achieves depar- 
tures from tracking, but no simple model has been exhib- 
ited in which the late-time behaviour necessarily approaches 
something close to a constant-density A term (Malquarti, 
Copeland & Liddle 2003). 

The other major problem with the quintessence ap- 
proach is that the models do not solve the scale problem: the 
potentials asymptote to zero, even though there is no known 
symmetry that requires this. This leads to consideration of 
the more radical many-universe mode of anthropic reason- 
ing. Here, one envisages making many copies of the universe, 
allowing the value of the vacuum density to vary between 
different versions. The simplest concrete way of generating 
this multiverse is via stochastic inflation with an additional 
scalar field that sets the effective value of A (e.g. Garriga & 
Vilenkin 2000). In this paper we will not need to be specific 
about the mechanism involved, although it is of course of 
the greatest interest. But first we must establish whether an 
ensemble approach makes sense. The logic is the same as 
in evolutionary biology, where we start with empirical evi- 
dence for selection out of a diversity of heredity, and only 
later move on to the search for the microscopic mechanism 
of genes and DNA that permits this diversity. 

Although most members of the hypothetical ensemble 
are presumed to have large vacuum densities comparable in 
magnitude to typical particle-physics scales, rare examples 
will have much smaller densities. Since large values of the 
vacuum density will inhibit structure formation, observers 
will tend to occur in models where the vacuum density falls 
in a small range about zero - thus potentially solving both 
the scale and why-now problems. The basic form of this ar- 
gument as an upper limit on the magnitude of A has existed 
for a number of years: see e.g. Linde (1984; 1987) or Bar- 
row & Tipler (1986). A significant further step was taken by 
Weinberg (1989), who used the anthropic argument to make 
the impressively bold prediction that A would indeed turn 
out to be non-zero at about the observed level. Weinberg's 
reasoning was taken up in more detail by Efstathiou (1995) 
and by Martel, Shapiro & Weinberg (1998). Efstathiou cal- 
culated the expected distribution for fly for a typical ob- 
server (a term whose meaning is discussed below); he found 
a result that peaked around flv — 0.9, in which the observed 
n„ = 0.75 would not be surprising. 

This is an impressive result, but it has two points at 
which further study is merited. Although he gave the gen- 
eral argument for the suppression of large observed values 
of A. Efstathiou then fixed the CMB temperature at its ob- 
served value in order to calculate a probability distribution 
for the observed value of f2„. In principle, it is possible that 
T — 2.725 K is not a typical value when all observers are 
considered, and we want to be sure that we have not brought 
the why-now problem in again by the back door. Garriga, 
Livio & Vilenkin (2000) have argued that anthropic selec- 
tion does indeed provide a general solution to the why-now 
problem, but it seems useful to make this argument explicit 
in terms of observables. One aim of this paper is therefore to 



revisit Efstathiou's calculation to calculate the joint distri- 
bution of A and T, to see how typical our observed location 
on this plane may be. We also take the chance to apply more 
accurate approximations for gravitational collapse of cosmic 
structures than the usual Press-Schechter (1974) approach, 
which underestimates the abundance of the most massive 
objects by a factor 10. 

A larger issue is the behaviour for negative A, which is 
examined in some detail. Efstathiou concentrated on models 
that were in the expanding phase, whereas a universe with 
negative A will eventually cease expanding and collapse to 
a big crunch. New accurate approximations are given for 
the growth of density fluctuations, for either sign of A. It 
turns out that structure formation in the collapsing phase is 
highly efficient, which presents something of a puzzle given 
the observed positivity of A. 

Finally, the same anthropic approach can in principle 
also be applied to the case of cosmological curvature. We 
indicate how this argument might have been applied to ex- 
plain an open universe in the period before evidence for A 
emerged. The paper concludes by discussing the testability 
of anthropic reasoning in the light of these results. 

2 ANTHROPIC SELECTION AND 
GRAVITATIONAL COLLAPSE 

2.1 Principles 

The first step in dealing with ensembles of universes is to 
decide what will be allowed to vary. In principle, nothing 
is guaranteed to be fixed, but it makes sense to start in as 
restrictive a way as possible. We therefore follow Efstathiou 
(1995): he assumed that all members of the ensemble are 
exactly spatially fiat, and that the key dimensionless ratios 
of cosmology are as we observe. These are (1) the photon to 
baryon number density ratio; (2) the dark matter to baryon 
density ratio; (3) the horizon-scale amplitude of density fluc- 
tuations, Sh- The only variation to be considered is in the 
value of the effective cosmological constant. 

Many authors have allowed a wider set of parameters to 
vary. For example, Garriga & Vilenkin (2006) consider joint 
variations in A and in 5h- The most radical view is that of 
the string-theory landscape, in which all of physics is free 
to vary (Susskind 2003), and aspects of this picture have 
been explored by e.g. Tegmark et al. (2006) and Graesser & 
Salem (2006). One approach to this question is an experi- 
mental one: if the simplest forms of anthropic variation can 
be ruled out, this might be taken as evidence in favour of 
the landscape picture. We therefore consider only the sim- 
plest picture here. This amounts in practice to considering 
the observed universe at some high temperature, so that any 
vacuum density of practical interest is negligible, and then 
considering the future evolution of copies of this universe in 
which the vacuum density is set to different levels. 

The most difficult and controversial question with any 
ensemble of universes is how to set a measure: what is the 
relative weighting of the different members of the ensemble? 
In the present context, we need to know the prior probabil- 
ity to assign to a given value of A, which will be modified by 
an 'observer-bias' factor, reflecting the relative difficulty of 
forming observers as a function of A. Various proposals for 
the measure in the case of eternal inflation have been consid- 
ered by Aguirre, Gratton & Johnson (2006). In the end, we 
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are persuaded by the original view of Weinberg (1989), who 
took the prior on A to be uniform over a small range around 
zero (see also Efstathiou 1995 and Weinberg 2000). If there 
is nothing special about a value zero, this seems a defensible 
assumption given that the range of interesting densities is 
tiny by comparison with particle-physics scales. 

Those different universes arc assumed to receive a 
weight according to the number of observers that exist in 
them, so that cosmologists asking questions here and now 
are treated as randomly selected from the totality of ob- 
servers over all universes and all times. The exact meaning 
of this intuitive idea of a 'typical' observer is not so eas- 
ily made precise, however, since observers are not standard 
objects (what relative weight do we assign to a modern hu- 
man, a caveman and a cat?). Some of these difficulties are 
reviewed comprehensively by Neal (2006). In the present in- 
stance, we can side-step many of these issues by exploiting 
the assumed similarity of the members of the ensemble in 
their non-vacuum physics. We do not need to predict the 
absolute number of observers, nor how they arc divided into 
different types of observer: it is sufficient to assume that a 
model with twice as many stars is twice as likely to be ex- 
perienced. Thus, we take the weighting of each member of 
the ensemble to be given by the fraction of the baryons that 
are incorporated into nonlinear structures. 

Weighting by collapse fraction is a common assump- 
tion in work of this sort, although not universally accepted. 
The most recent dissenters are Staxkman & Trotta (2006), 
who calculate the maximum number of observations a given 
observer can ever make, which introduces a further depen- 
dence on A. There is certainly an ambiguity worth debating 
here: given two universes that make the same number of ob- 
servers, but where one set live much longer than the other, 
should we give an equal weight to each, or weight them in 
proportion to lifetime? Starkman & Trotta effectively treat 
equal time intervals as equally probable, but the instants in 
the life of a given individual are not independent: any of us 
will always answer "yes" to the question, "is A surprisingly 
low?" , and it should not matter how many times we are able 
to ask. Therefore, we prefer to weight by the numbers of ob- 
servers produced, and it seems reasonable to tie this to the 
number of stars that are formed. 

As discussed earlier, we need a prior for the different 
values of the vacuum density: following Efstathiou (1995) 
and Weinberg (1989; 2000), this is taken as being uniform 
over a range around zero. We therefore weight universes with 
an overall posterior probability 

dP{pv) oc fc dpv, (2) 

where fc is the collapse fraction: the proportion of mass in 
the universe that has become incorporated into sufficiently 
large nonlinear objects. 

In more detail, we will be interested in the weight at- 
taching to different times in each universe. This is in part 
a simple issue, since we can calculate what fraction of the 
mass undergoes collapse in any given time interval, so that 

dP{pv,t) oc ^ dt dpv (3) 

But this gives the time distribution for the formation of sites 
at which life might subsequently emerge. The more serious 
challenge lies in predicting the history of observers follow- 



ing a formation event: how many observers will eventually 
result, and what will their distribution in time be? Rather 
than making arbitrary assumptions, we can avoid the worst 
uncertainties by turning the problem backwards. We can 
calculate the distribution of times at which stars form in 
the universe, and we know when the star with which we are 
associated was formed: 4.6 Gyr ago. That time corresponds 
to a redshift 0.457, at which point the cosmological param- 
eters were T = 3.97 K and n„ = 0.49. We can therefore 
concentrate on the more concrete question of whether the 
sun formed at a typical point in comparison to all stars in 
the multiverse. This presumes that the subsequent history 
of life ajound each star has no dependence on A, but further 
biological assumptions are not required. 

2.2 The collapse fraction 

For this calculation, we need to be able to predict the frac- 
tion of baryons in the universe that are processed into stars. 
This is not something that can presently be calculated with- 
out some guidance from observation. The first galaxy-sized 
systems to collapse are of low mass and high density, and 
clearly will form some stars, but it is expected that feedback 
from this initial activity will quickly regulate the star forma- 
tion in these objects. Thus in practice most star formation 
is expected to occur in the largest galaxies, and we follow 
Efstathiou (1995) in treating these as being defined by a sin- 
gle mass scale. However, as shown below, it is now possible 
to deduce this scale empirically, rather than appealing to a 
priori cooling arguments such as Rees & Ostriker (1977). 

We develop below the formulae needed to calculate ac- 
curately the fraction of the mass that has collapsed into 
objects of a given mass scale or larger. As usual, the mass 
scale is defined by the mass in a homogeneous universe con- 
tained within a sphere of radius R. The fractional density 
fluctuations smoothed with such a spherical filter have an 
rms value a{R), and the rareness of objects of a given mass 
is quantified by defining 

v = 5rJcj{R), (4) 

where 5c is a density threshold of order unity. Since a 
changes with time, we need to specify an era in order to 
associate p with a given mass. It is convenient to make the 
arbitrary choice of T = 1000 K as a reference era (matter 
dominates over radiation and over any vacuum density of 
interest). Fig. 1 then shows that the existing data on the 
evolving comoving stellar density can be well described us- 
ing a single-scale model with 

v{T = 1000) = 250, (5) 

adopting 5c = 1.686 as justified below. 

This plot does not need to make any assumptions about 
the matter power spectrum, but if we adopt a canonical 
model (O^ = 0.25, fif, = 0.04, h = 0.73, Us = 0.95, as = 
0.8), this corresponds to an effective galaxy mass of 

Mg = 1.9 X lO^^M©. (6) 

The small-scale CDM power spectrum is rather flat: 

dlni^/dlnM~ 0.145, (7) 

so the results of this calculation should be relatively robust. 
Loeb (2006) is correct to point out that, in principle, some 
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redshift z 

Figure 1. Time dependence of star formation predicted in a sim- 
ple collapse model. The total stellar density produced by a given 
epoch is assumed to scale with the total collapse fraction associ- 
ated with a single mass scale. The density-fluctuation parameter 
i^(T = 1000) is varied by up to 20 per cent either side of its canon- 
ical value u = 250. This yields a good match to the data on the 
empirical redshift dependence of the total stellar mass density, 
taken from Merloni, Rudnick & Di Matteo (2004). 

stars could form in dwarf galaxies, but changing the critical 

mass scale by a power of ten has little impact on the results. 

Even so, it is clear that this simple model can be chal- 
lenged, since the typical galaxy mass arises in a complex way, 
where feedback from supernovae and AGN heats gas and 
prevents star formation being very efficient (the observed 
cumulative efficiency today is flt/ilb — 5% — 9%, depending 
on assumptions about the IMF: Cole et al. 2001). It is cer- 
tainly possible that the operation of such effects in objects 
of a given mass could depend on density and thus on era - 
so that the overall efficiency of star formation could have a 
further complicated dependence on A. We shall ignore this 
point here, but it would clearly be of interest to investigate 
the issue elsewhere using detailed galaxy formation models. 



2.3 The cosmological mass function 

Given a filtering scale corresponding to a typical galaxy, the 
linear density contrast, cr{R) can be calculated. According 

to the Prcss-Schechter (1974) approximation, the collapse 
fraction (i.e. proportion of mass contained in objects of the 
given mass scale or larger) is then 

fc = erfc {u/V2) ■ V = 5c/ cr. (8) 

The critical density contrast is 5c — 1.686 in an Einstcin-de 
Sitter model. Efstathiou (1995) gives an approximate scal- 



ing as (1 — n^) 



but this is a theoretical expectation 



based on the spherical collapse model. In detailed studies 
of numerical simulations, Jenkins ct al. (2001) found that 
5c could be treated as constant in matching theory to the 
simulated mass functions. Jenkins et al. followed Sheth & 
Tormen (1999), who established that the Press-Schechter 
form for the mass function was significantly in error, with 
too many objects at the peak of the mass function by about 



a factor 1.5, and too few at the highest masses by a power of 
10. More recently, Warren et al. (2006) have shown that the 
Jenkins et al. fitting formula still contains errors of order 
10%, and they proposed a replacement. An unsatisfactory 
feature of their fit is that it predicts collapse fractions in 
excess of unity for small masses, whereas one would nor- 
mally prefer to assume that /c — > 1 in the limit of very 
small masses. It is simple to cure this by finding an ana- 
lytic formula for fc; this can then be differentiated in order 
to find the mass function. The following expression matches 
the Warren et al. fitting formula to a maximum error of 
about 1% over the whole range where data exist: 

/c = (l+oi^'')-'exp(-c;/'), (9) 

where (o,6,c) = (1.529,0.704,0.412). 



2.4 Growth of density perturbations 

For this exercise, we also require the linear growth func- 
tion for density perturbations, which can be expressed as 
a numerical integral (Heath 1977). It is convenient to have 
an accurate numerical approximation, and the following ex- 
pressions are good to a maximum error of 0.1%. The cases 
of positive and negative A are somewhat distinct. For the 
positive case, 

5{a) ~ x(l - x^-^y-^^ + 1.437 (l - (1 - x^^^) , (10) 

where x denotes fli, (a)^'''^, and we choose the a = 1 point to 
correspond to equal density in matter and vacuum: 

-3\-l 



Q.v{a) = (1 -(- a 



(11) 



so that 5{a) ~ a for small a. For a starting point where 
Q.V (o) is small, the total amount of growth is 



5oo/5iT, 



1.437/ainit ^ 1.437/[n^(oinit)]'''^ 



(12) 



Since f2„(ainit) = Pi,/5.375 x 10^ moV* for our choice of 
T = 1000 K as a normalization point, at which the galcixy- 
scale fluctuation is = 250, this immediately allows the 
asymptotic value of v to be deduced. 

For the negative density case, we need time as a coor- 
dinate, since the scale factor is not monotonic: 



a{t) = [sin(3t/2)]2/^ , 



(13) 



where hero we choose units such that a = 1 at the point 
of maximum expansion, and time is measured in units of 
(87rG|p„|/3)~^^^, so that Friedmann's equation is (a/a)^ = 
— 1 and 0^(o) = (1 — a~^)~^. Here, the approximation 
for the growth function is 



m = 



(14) 



(1 + 0.37(t/icoii)2-i*) (1 - (t/icoii)2) ■ 

Again, the normalization is that 5{a) ~ a for small a. Note 
that the fluctuations diverge at the collapse time (tcoii = 
27r/3) as l/(icoii —t): this corresponds to the decaying mode 
in the expanding phase. 

3 THE EXPECTED SIGN OF A 

The collapsing phase turns out to be important when we 

consider structure formation. A substantial negative A lim- 
its the amount of growth that can occur before the universe 
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Figure 2. The collapse fraction as a function of the vacuum 
density, which is assumed to give the relative weighting of different 
models. The dashed line for negative density corresponds to the 
expanding phase only, whereas the solid lines for negative density 
include the recoUapse phase, up to maximum temperatures of 
10 K, 20 K, 30 K. The observed value of the vacuum density is 
33meV4. 

ceases to expand, but the total amount of growth after this is 
Umitod only by how close to the big crunch we are prepared 
to venture. Eventually, /c tends to unity in all such recoUaps- 
ing models. Since the uniform prior extends arbitrarily far 
towards more negative values of A, this apparently implies 
that all the weight should be given to models with A < 0. 

However, structures that form very close to the final 
singularity are not of interest for the anthropic calculation: 
there is little time remaining for life to develop, and in any 
case the CMB will have heated up to the point where it 
interferes with life - or indeed perhaps even with the forma- 
tion of stars and planets themselves. It is simplest to express 
this cutoff in the recoUapsing phase in terms of a mgiximum 
temperature that we are willing to consider, although this 
can be directly translated to a limit on the time remaining 
before the big crunch. Since the recoUapsing phase is the 
time-reversed version of the expansion, the time remaining 
from temperature T until the big crunch is just what would 
have elapsed from the big bang until this temperature. Nor- 
mally, the matter-dominated approximation will apply, so 



tiT) 



3Ho 



-1/2 



(rp X -3/2 



(15) 



We know from observations that star formation in galaxies 
can proceed actively at redshift ~ 7, so Tmax > 10 K on 
these grounds. This would leave only a few Gyr after for- 
mation for life to evolve, so presumably Tmax should not 
be much larger than this, and could well be smaller. This 
is not so much a biological argument as one based on stel- 
lar lifetimes. Highly negative values of A would cause the 
universe to turn round without ever cooling below this criti- 
cal temperature, so these models may be excluded from the 
point of view of generation of observers, even though they 
form galaxy-scale structures efficiently. This is completely 
distinct from the situation at highly positive A, where the 
problem is the failure to create structure at any time. 



Rather than singling out a particular value of Tmax, 
we may as well perform the calculation allowing it to take a 
range of values. This is illustrated in Fig. 2, which shows the 
relative weight to be given to models as a function of p„ . We 
can integrate this distribution to obtain Fig. 3, which shows 
how the probability of observing a negative A varies with the 
maximum temperature for structure formation that we are 
willing to tolerate. The probability that observers inhabit 
a universe with A < is about 50% if Tmax = 8.5 K, or as 
little as 1% if Tmax ~ 1.5 K. Conversely, only about 6% of ob- 
servers will experience a positive A if Tmax = 30 K; but this is 
probably too tolerant of late- forming observers. From Fig. 2, 
it is furthermore clear that the bulk of any weight in favour 
of negative A lies with the recoUapsing phase: the probabil- 
ity of inhabiting the expanding phase is approximately the 
same as the probability that A is positive. According to the 
anthropic framework, it is therefore far from inevitable that 
we have ended up in a positive-vacuum expanding universe, 
although neither is it particularly unusual. 

We can now make a first attempt to assess how well this 
anthropic prediction matches reality. Since A = is a special 
point, it is reasonable to consider the frequentist question: 
'what is the probability that lies within (2.39 meV)"' of 
the origin?'. The answer is also plotted in Fig. 3: about 10% 
if we ignore negative A, peaking at 20% for Tmax — 4K, and 
declining for larger values. The consistency of observations 
with anthropic prediction therefore depends somewhat on 
the recoUapsing phase. In what follows, we will concentrate 
on the expanding case of positive A, but it should be borne 
in mind that any anthropic probabilities we deduce subse- 
quently should be reduced slightly to account for the fact 
that observers also have a non-negligible chance of being 
found in a recoUapsing model. 



4 THE JOINT DISTRIBUTION OF A AND T 

The above calculation addresses the scale problem, but only 
considers the total amount of structure formation, not when 
it occurs. This problem was considered by Garriga, Livio & 
Vilenkin (2000), who calculated the posterior distribution 
on the plane (tAjio), where tA is the time of A-domination 
and tc the time at which a typical galaxy was formed. They 
used a Press-Schechter approach to show that the proba- 
bility density on this plane peaks around to ~ ^a, so that 
the why-now problem is solved. We need something similar, 
but we want to see explicitly how the Efstathiou analysis 
is affected by a given assumed temperature; it is therefore 
necessary to cast the distribution in terms of the observables 
n„ and T. We have the differential probability distribution 
dP oc dpv dfc, and we want to change variables to Civ and T. 
Rather than doing this in a single step by working out the 
Jacobian of the transformation, we can first note that the 
collapse fraction is just a function of T for given so that 



dP oc dpv 



dT. 



(16) 



If we now change from {pv,T) to (n„,T), the Jacobian is 
diagonal, and 



dP oc 



dQ,v 



T dT 



dQv dT. 



(17) 
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Figure 3. Integrating under the distribution of Fig. 2, we can 
deduce the probability of inhabiting a universe with A < 0, as a 
function of the assumed maximum temperature for galaxy forma- 
tion (solid line). The dashed line shows, under the same assump- 
tion, the probability that the absolute value of A would lie within 
its observed value of (2.39 me V)*. For a maximum temperature 
of 8.5 K, positive and negative values of A are equally probable, 
but negative A is disfavoured by more stringent limits on tem- 
perature: for a maximum temperature of 1.5 K, only about 1% of 
observers will inhabit a universe with A < 0. 




Figure 4. The probability distribution of Qy at observed temper- 
atures of 1 K, 2 K, 4 K, 8 K, with higher temperatures pushing the 
distribution to lower values of Qy. This plot should be contrasted 
with Fig. 2 of Efstathiou (1995). 



To evaluate the first factor on the rhs, note that the vacuum 
density parameter at some early time to is 



Qv + {i — iiv)aQ (1 — fi^)oo' 



(18) 



where fly is the density parameter at the later time of in- 
terest, and ao = T/Tq oo gives the latter approximation. 



In this limit, r2„o -C 1; this is therefore proportional to p„, 
giving 



(19) 



This gives higher weight to Q.^ close to 1, since this is an 
attractor for the evolution of Q,v{t) if A > 0. The focusing 
towards Q.v = 1 increases as T falls, and this is reflected in 
the factor. As a result, dP/dT with n„ fixed at zero difi^ers 
in this framework from what we would normally calculate 
for an Einstcin-do Sitter universe. 

Efstathiou (1995) did not need to consider the T'^ fac- 
tor, since he held the temperature constant at its observed 
value. This is appropriate if we want to use the anthropic 
framework in a Bayesian sense, to make the best predic- 
tion of the current value of r2„ given what else we know. 
For similar reasons, Efstathiou held constant the observed 
large-angle CMB fluctuations, whereas these vary with tem- 
perature (i.e. with time of observation). Here, we axe inter- 
ested in the broader question of whether the conditions we 
observe are close to those experienced by a typical observer. 
To show that this distinction matters, consider Fig. 4. This 
plots the posterior probability of Q.v for various choices of the 
observed temperature, and shows that the result is sensitive 
to temperature. For T ~ 8 K, the distribution peaks near 
the observed Q.v ~ 0.75, but for lower T the distribution is 
dominated by the spike in the prior at n„ = 1. 

We now return to the full joint distribution for f2„ and 
T, which we had in the form 



dP oc 



dO,v dT. 



The remaining partial derivative is 



dfc 
dT 



dlnf d\nT' 



(20) 



(21) 



For convenience, we can use the Peebles (1980) approximar 
tion for the logarithmic growth rate: 



91ni^ 
d\nT 



so that overall 

dfc 



dP(x 



dim 



T^(l 



dQv dT. 



(22) 



(23) 



This joint probability distribution is show in Fig. 5, con- 
verting to log T for convenience. The 'observed' universe of 
{T,fly) = (3.97,0.49) is plotted as a point. It is clear that 
the point does not lie in a particularly unusual position in 
this plane. If wc draw contours of constant probability den- 
sity, the point lies well within the 68% contour. These con- 
tours are not unambiguous, as they depend on the measure 
adopted on the Qv — T plane. However, if wc inspect the 
marginalized distributions for $1„ and T, shown in Figs 6 
& 7, we see that the observed conditions are close to the 
50% point in each quantity. In short, the anthropic calcula- 
tion suggests that we are indeed extremely typical observers, 
both in terms of the vacuum density we see, and when we 
see it. 
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T / K 

Figure 5. Contours of probability density on tlic log(T) — 
plane. The three contours shown enclose 68%, 95% and 99% of 
the probability. The solid point shows the conditions at the epoch 
of formation of the sun. 




Figure 6. The integral probability distribution for f2„, marginal- 
ized over T. 



5 ANTHROPIC ARGUMENTS APPLIED TO 
CURVATURE 

So far, we have assumed that the vacuum is indistinguish- 
able from a cosmological constant, with equation of state 
w — P/pc = — 1. Wc know that this is a good approxima- 
tion for our universe, and it would in any case go beyond 
the scope of this paper to consider ensembles in which more 
than one parameter varies. However, it is worth paying some 
attention to the special case w = —1/3. With the critical ex- 
ception of the distance-redshift relation, a flat model with 
this vacuum equation of state is indistinguishable from a 
model with non-zero spatial curvature: the Priedmann equa- 
tion and the growth equation for density perturbations are 
identical. We can then use a modified version of the above 




T / K 

Figure 7. The integral probability distribution for T, marginal- 
ized over ilv 

approach to show what happens if we confront anthropic 
reasoning with the curvature of the universe. 

To some extent, curvature presents a parallel set of 
problems to the vacuum. There is a scale problem, in the 
sense that natural initial conditions might be thought to 
have a total |f2 — 1| of order unity, which would lead to a 
universe dominated by curvature long before today. There 
could also be a wliy-now problem, if the present curvature 
was non-zero at the level of |f2 — 1| ~ 0.01, which cannot 
currently be excluded. It is commonly assumed that infla- 
tion solves the curvature scale problem and also predicts 
that there is no why-now problem, but it is interesting to see 
how the anthropic apparatus copes with this case. This issue 
has been examined previously (e.g. by Garriga, Tanaka & 
Vilenkin 1999), but it is of some interest to see how our spe- 
cific approach works out in this case. Rather than break our 
rule of allowing only one parameter to vary in the ensemble, 
we take a historical approach and imagine how anthropic 
arguments might have been applied to curvature decades 
ago, when many cosmologists were convinced that A = 
(it is perhaps surprising that anthropic ideas received little 
emphasis at this time). 

Unlike the vacuum density, curvature lacks an obvi- 
ous time-independent absolute scale. At any given era, one 
can define Qk = 1 — f2m, so it will be convenient to con- 
sider r2fc(r — 1000) as our parameter. This evolves as 
f2fe(a) = ^k/{^k +f2m/a)- Interesting values of this number 
at the reference T = 1000 K will be small, so it is tempt- 
ing to follow our previous procedure and assign a uniform 
prior around zero, and weight models by their asymptotic 
collapse factor. In the days before inflation, this might have 
been a defensible expression of ignorance: the essence of the 
flatness problem is that order unity positive or negative cur- 
vature in the initial conditions seems more natural than the 
tiny amount necessary to yield an almost flat universe today. 
But in modern models where 'pocket' universes are formed 
by tunnelling, the result is an open universe, so priors on cur- 
vature might well have a discontinuity at zero (e.g. Preivo- 
gel et al. 2006). The idea of a uniform prior for curvature 
is therefore less well founded than it is for A. Nevertheless, 
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it is of some interest to carry out the exercise of adopting a 
uniform prior and seeing where it leads. 

For models with A = 0, the perturbation growth as a 
function of a is analytic: 



6{a) = +1 + — - vTT^sinh-^v^) (0^ > 0) 



-1 + 



a3/ 
3 

,3/2 



(\/a — \/l — asin ^\/aj (fife < 0); 



(24) 



in the collapsing phase for negative Qk, sin~^y^ is replaced 
by sin~^v^ — tt. These expressions are normalized so that 

S ~ 2a/5 for a <^ 1. The convention for the scale factor 
assumes that a can be written in terms of conformal time, 
77 as 

o = (1 — cosr))/2 (closed) 
a = (coshr; — l)/2 (open), 



(25) 



so that a = 1 at ma^ximum expansion in the closed model. 
The relation to density parameters in these units is 



f^fc(a) = a/ {a — 1) (closed) 
fifc (a) = a/(a + 1) (open). 



(26) 



We can now repeat the exercise of Fig. 2 for the case 
of curved universes with A = 0, and the results are shown 
in Fig. 8. Generally, the anthropic weighting as a function 
of curvature looks similar to the weighting as a function 
of A, but with some important differences. The asymme- 
try in favour of recoUapsing models is not so extreme: the 
probability of experiencing 0^ > (i.e. a negatively curved 
open universe) is 41% for Tmax = 10 K, falling to 29% for 
Tmax = 30 K. 

What about the magnitude of curvature? A decade ago, 
open models were seriously under consideration, and some 
would have argued for Q.^ — 0.7, so that ilk{T = 1000) ~ 
0.006. From Fig. 8, we see that the typical curvature pre- 
dicted for such A-free universes was ^^(r = 1000) ~ 0.01, 
so an anthropic approach to explaining the density param- 
eter in matter-only models would have yielded sensible an- 
swers. In order to reject an anthropic explanation for an 
open universe at the 1% level, it would have been necessary 
to have a limit of \QkiT = 1000) | < 10"" (depending slightly 
on limiting temperature), corresponding to \Qk\ 5 0.035 at 
T = 2.725. In fact, we barely know that the universe is flat to 
this precision even today: the present limit is approximately 
\Qk\ < 0.02, according to Spergel at al. (2006). Therefore, an 
anthropic approach to curvature would have been perfectly 
consistent with 1990s data. 

Today, the issue of curvature would be approached in 
the context of inflation, where a sufficient number of e- 
foldings, iV, of the expansion will lead to values of the 
present curvature that are unmeasurably small. The issue 
of interest is therefore the prior to be placed on A'^ (see e.g. 
Freivogel et al. 2006) . Whatever the result of such a calcula- 
tion, however, a sufficiently small upper limit on curvature 
would allow the anthropic argument to be rejected. Simi- 
larly, if we had no detection of A, our earlier results show 
that a sufBcicntly strong upper limit would reject the an- 
thropic approach, leading us to require a physical mecha- 
nism that forces A = 0. Anthropic reasoning is thus testable 
and could point to new physics. But this is not the situation 
we face: we have an actual detection of A, rather than an 
ever-retreating upper limit, and no a priori theory predicts 



the observed number. An explanation in terms of anthropic 

selection from an ensemble matches what we see, and so far 
there is no credible alternative. 



6 CONCLUSIONS 

This paper has revisited the anthropic calculation of Efs- 
tathiou (1995) in more detail, dropping the assumption that 
the CMB temperature is fixed, and considering also the for- 
mation of structure during the collapsing phase of a model 
with negative vacuum density. We adopt Weinberg's (1989) 
assumption that the prior on A is uniform in a small range 
around zero, so that models are weighted simply by their col- 
lapse function. We have attempted to avoid issues such as 
the lifespan of civilizations by using the formation era of the 
sun as the data to explain: so long as observers eventually 
form with some fixed mean number per star, we can ignore 
how they are subsequently distributed in time. For universes 
that are expanding, the conclusion is that the sun formed 
at completely typical values of n„ and the CMB tempera/- 
ture. There is therefore no observational reason to challenge 
the idea of a uniform prior on A, nor to require variation of 
other parameters within the ensemble although there re- 
mains theoretical motivation for considering more complex 
alternatives. 

This approach can certainly be questioned in the case 
of negative vacuum densities, where the age of the universe 
is finite. One might bo tempted to argue that all the weight 
should go to positive-density universes, since civilizations 
can then potentially last forever. However, this is too opti- 
mistic, since the event horizon in models that are asymptot- 
ically de Sitter limits the resources that are available. Once 
most of the stars in a given model have formed, it is reason- 
able to expect that observers will find existence progressively 
more difficult after a further few billion years as the existing 
stars die out. In any case, we know that we are members of a 
civilization that has not yet outlived its star; the question of 
whether some races might live for a trillion years can there- 
fore be ignored. The reasoning here is the same as in the 
"God's coin toss" thought experiment discussed by Olum 
(2002). This experiment imagines that according to the toss 
of a coin either 10 (heads) or 1000 (tails) people are cre- 
ated and given a number. If you have no knowledge of your 
number, your odds for the result of the coin toss should be 
100:1 in favour of tails; but if your number is < 10, your 
odds should be equal. 

From this point of view, there thus seems no reason not 
to give appropriate weight to the flourishing of observers 
in recoUapsing models. In contrast to models with positive 
vacuum density, gravitational collapse is always perfectly 
efficient in models whore A < 0. But this docs not imply 
that anthropic selection must favour this case: temperatures 
should not too be too high at putative formation, and there 
should remain at least a few billion years before further col- 
lapse renders the CMB hot enough to be an environmental 
hazard. These criteria motivate a maximum temperature of 
order 10 K, for which only a minority (~ 10%) of observers 
should find themselves in recoUapsing models. This was not 
our fate; but it is interesting to speculate how observational 
cosmology might have developed in such a case. 
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Figure 8. The collapse fraction as a function of the curvature 
in models with A = 0, which is assumed to give the relative 
anthropic weighting of different models. Curvature is specified as 
n^. = 1 — at the era when T = 1000 K. The dashed line for 
negative corresponds to the expanding phase only, whereas 
the solid lines for negative Oj. include the recoUapse phase, up to 
maximum temperatures of 10 K, 20 K, 30 K. 
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